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Abstract 

Within the framework of the operator product expansion (OPE) and the 
renormahzation group equation (RGE), we show that the temperature and 
chemical potential dependence of the zeroth moment of a spectral function 
(SF) is completely determined by the one-loop structure in an asymptotically 
free theory, and in particular in QCD. Logarithmic corrections are found to 
play an essential role in the derivation. This exact result constrains the shape 
of SF's, and implies striking effects near phase transitions. Phenomenological 
parameterizations of the SF, often used in applications such as the analysis 
of lattice QCD data or QCD sum rule calculations at finite temperature and 
baryon density must satisfy these constraints. We also explicitly illustrate in 
detail the exact sum rule in the Gross-Neveu model. 
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I. INTRODUCTION 



Hadronic properties at finite temperature and baryon density are of great importance in 
tlie pfienomenology of lieavy ions collisions, star interior and the early universe. Moreover, 
the theoretical expectation of transitions to a chirally symmetric phase and, perhaps, to a 
quark-gluon plasma phase contributes to the interest in studying the effect of matter and 
temperature on the quantum chromodynamics (QCD) vacuum. 

Our present understanding of QCD at finite temperature (T) and baryon density (or 
chemical potential /i) is mainly limited in the Euclidean realm, due to the lack of non- 
perturbative and systematic calculating tools directly in the Minkowski space. Typical 
methods, with QCD Lagrangian as the starting point, are the OPE and lattice simula- 
tions. Because these two formulations are intrinsically Euclidean, only static quantities are 
conveniently studied. In order to gain dynamical informations, which are more accessible 
experimentally, the analytic structure implemented through dispersion relations often have 
to be invoked within the theory of linear response. 

In principle, dispersion relations allow the determination of the spectral function (SF), 
which carries all the real-time information, from the corresponding correlator in the Eu- 
clidean region. In practice, realistic calculations, e.g. OPE or lattice simulations, yield only 
partial information on the correlators, making impossible a direct inversion of the dispersion 
relation. 

Therefore, the standard approach assumes a phenomenological motivated functional form 
with several parameters for the SF, and uses the information from the approximate corre- 
lator, mediated by the dispersion relation, only to determine the value of parameters by 
a fit. This approach has been quite successful at zero temperature and density, thanks to 
the fortuitous situation that we roughly know how to parameterize the SF's in many cases. 
Two important examples are the QCD sum rules pioneered by Shifman, Vainshtein and 
Zakharov lll,^! and the analysis of lattice QCD data [Q. So far, standard parameterizations 
have included poles plus perturbative continuum [Q. 

The success of such approaches heavily rests on our good understanding of the qualitative 
behavior of SF's at zero (T, /i). We can find other such favorable examples in the low-T 
regime , where the shape of the SF has the same qualitatively features of the zero (T, /x) 
case, or even in the high-T regime for simple models for which the functional form of 
the SF is known. 

The QCD sum rules approach has been extended also to systems at finite tempera- 
ture The lack of experimental data, and of reliable nonperturbative calculations has 
prompted people to use the same kind of parameterizations that have worked so well at 
zero temperature with, at most, perturbative corrections. We believe that physical results 
at finite (T, /i) can be strongly biased by this assumption. In fact, recent interpretations of 
lattice simulation data P-[TT[] appear to indicate the existence of such problems. 

The purpose of this work is to derive exact sum rules that constrain the variation of 
SF's with (T, /i). In addition, we apply these sum rules to the chiral phase transition, and 
demonstrate that SF's in some channels are drastically modified compared to both their 
zero (T, fi) and their perturbative shapes. This result confirm our worries about non-trivial 
effect of finite T or baryon density on the shape of the SF. 

Our derivation of these exact sum rules, based on the OPE and the RGE, has a closer 
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relation with the derivation of sum rules for deep inelastic scatterings than with the QCD 
sum rule approach of SVZ In fact, we establish relationships between moments of the 
SF and corresponding condensates as functions of (T, /i), without assuming any functional 
form of the SF. 

In the derivation process, we find that the logarithmic corrections are essential to estab- 
lish the exact sum rules. In contrast, the QCD logarithmic corrections are only marginally 
relevant in the finite energy sum rules, and hence are rarely discussed in the literature. To 
properly take into account the logarithmic corrections, a repeated partial integration method 
is used to match the relevant asymptotic expansions. 

Since no further assumptions other than the validity of the OPE and RGE are involved 
in the derivation, our sum rules are very general and can be applied anywhere in the (T, /i)- 
plane, even near or at the boundary of a phase transition. 

The paper is organized as follows. In section II we present the general derivation of exact 
sum rules in asymptotically free theories. The matching of asymptotic behaviors of the 
correlator and the dispersion integral, including explicitly their logarithmic corrections, are 
carefully discussed. In section III we illustrate each single step of the derivation in a soluble 
model, the Gross-Neveu model in the large- limit. In this model we can calculate exactly 
all the relevant quantities (spectral functions, Wilson coefficients, condensates, anomalous 
dimensions and correlators in space-like region, etc.), and, therefore, give a concrete example 
of how our method works. The application of our method to the derivation of exact sum 
rules for the mesonic channels in QCD is presented in section IV. In the same section, we 
also discuss the phenomenological consequences of the exact sum rules near chiral restoration 
phase transitions. Finally, in section V we summarize our work, draw some conclusions, and 
discuss possible future directions. 

II. GENERAL DERIVATION 

We start this section with a short review of the linear response theory, the OPE and the 
RGE. Next we introduce a convenient subtraction in the dispersion relation for studying the 
dependence of SF's on T and /i. 

Then, we present in detail a crucial part of our method: how to match the asymptotic 
OPE expansion with a corresponding asymptotic expansion of the SF and its dispersion 
integral. This approach is necessary for properly taking into account the logarithmic cor- 
rections, and studying the convergence properties of the relevant moments of the SF. More 
naive approaches not only yield, in general, incorrect sum rules, but might also fail to rec- 
ognize that a given sum rule does not exist in the first place, since the integral of the SF 
involved is divergent. 

Finally, a comparison of the two asymptotic expansions leads to the desired exact sum 
rules. We end this section with some general comments on the derivation and meaning of 
these sum rules. 
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A. Linear Response 

The real-time linear response to an external source, S{x), coupled to a renormalized 
current J{x) in the form of AC = S{x)J{x) is given by the retarded correlator: 

Kix;T,f^)^9ixo){[Jix),J{0)])T,,, (1) 

where the average is on the grand canonical ensemble specified by temperature T and chemi- 
cal potential /i. Disregarding possible subtraction terms, which are (T, /i)-independent since 
they are related to short distance properties of the theory, we can write the following dis- 
persion relation for the frequency dependence of the retarded correlation function: 

K{u},k;T,fi) = du — — — . (2) 

JO — [uj + ley 

For convenience, we discuss only the uniform limit {k = 0) and, from now on, drop the 
momentum label. We wish to emphasize that, even if the derivation of the sum rules requires 
the analytic continuation of Eq. into the Euclidean region {uj —>■ iQ), the dispersion 
integral is, nevertheless, determined by the singularities of p{u; T, /z) in the Minkowski region. 
One should keep clearly in mind the distinction between these real-time singularities with 
the ones related to screening phenomena. 



B. The OPE and the RGE 



In asymptotically free theories, the correlator K{uj = iQ; T, fi) is calculable in terms 
of the OPE [|13[ asymptotic series in the limit of oo. For a given renormalization 



prescription, the asymptotic expansion reads 

K{tQ, k; T, fi) ~ KoitQ, «:) + ^ C„(Q', m\f,) , g\t,) , «:) ([0„]«)t,^ , (3) 

n 

where m^K), q^k), the [0„]k's and the C„'s are, respectively, the current quark mass, the 
coupling constant, the renormalized composite operators, and their corresponding Wilson 
coefficients at the subtraction scale k. Here and in the following, we use the symbol ~ 
to denote asymptotic equality. We explicitly singled out the perturbative term Ko^iQ,/^), 
which corresponds to the identity operator and, being independent of T and fi, is calculated 
at {T,fi) = 0. This term is the only one that is not suppressed by powers of l/Q"^ and 
can depend on an overall subtraction scale, which we identify, for simplicity, with the same 
subtraction scale k used in the renormalization of the operators. It is important to note 
that the information of the ensemble average is encoded in the matrix elements of composite 
operators, while the Wilson coefficients, including Kq, are independent of T and fi. Although 
the matrix elements {[On]K)T,fj. cannot be determined perturbatively, the Q^-dependence of 
the Wilson coefficient C„ is controlled by the following renormalization group equation 



K^ + 2Tj-T„]CniQ^m\K),g\K),K) = 0, (4a) 
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where 

d d ^ d ^ d , 
an OK og om 

and Fj, r„ and Tm are the anomalous dimensions for the current J, the operator On 
and the current quark mass m respectively. For the purpose of illustration, we are only 
considering operators that do not mix under renormalization, but the mixing will be properly 
taken into account when necessary. The standard approach |]T^ to the renormalization 
group equation is the introduction of a running coupling g^iQ) and a running mass m{Q). 
In asymptotically free theories g'^{Q) vanishes logarithmically at large Q^. It is therefore 
meaningful to consider, in this limit, a perturbative expansion of the renormalization group 
equation functions 

r, = -7,^2 ^ ^^gA^ ^ and /5 = -bg' + 0{g') , (5) 

where i = J, n, m, while 7j, b are pure numbers determined by a one-loop calculation. 
Within this perturbative context, Eq. can be solved |TB| : 



C„(g^m(K),,^(«:),.)=(i^)^ (^) " xc„(^,,^(Q)), (6) 

where dn is the canonical dimension of the composite operator On minus the dimension of K 
in units of a mass, and the Cnijn? {Q) / Q"^ , fl'^(Q))'s are calculable perturbatively. Corrections 
to this result coming from nonleading powers of g^ in the expansions of the renormalization 
group equation functions, see Eq. (|]), are analytic in g^{Q) = 1/(6 In Q^) and hence can be 
absorbed into c„. Therefore, the leading term in the large-Q^ limit is totally characterized 
by the one-loop structure of the theory. 



C. Subtraction and (T, /i)-dependence 

To study the temperature and chemical potential dependence of i^T , we only need to 
consider the difference AK{iQ): 

AKm = K{zQ, K- T, fi) - k{iQ, K- r, fi') = du^ , (7) 

where Ap{u) = p{u;T,fi) — p{u; T' , p') . This subtraction is crucial to remove Ko{iQ,n), 
which contains terms not suppressed by powers of and is explicitly dependent on the 
renormalization point. The OPE asymptotic expansion of AK{iQ) is then 

^'<m ^ X .„,!^„^,«),A(|o.u (8) 

where A([0„]k) denotes the difference of the expectation values of [O^]^ in the ensembles 
specified by (T, /i) and (T', /i'), respectively. Since the current quark mass only runs logarith- 
mically in asymptotically free theories, we can safely ignore the corrections to the c^'s due 
to their dependence on m?{Q)/Q'^ when in deep Euclidean limit. Notice that the subtrac- 
tion has made AK{iQ) independent of the renormalization point k, i.e. AK{iQ) satisfies a 
homogeneous RGE. 
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D. Asymptotic Expansion 



Since Cn{m? {Q) / Q"^ , Q^iQ)) is perturbative and hence can be expanded in power of g^{Q) 
(and denoting the first non- vanishing power as z^^"-*), the left-hand side of Eq. can be 
expressed as a double (in and g^{Q)) asymptotic expansion of the form: 

ri,,i/=0 ^ 

where the exponents ?7„ = + {2'jj — 7^)/ (26) and the Q^-independent coefficients c^''-'(k) 
are again known perturbatively. It is important to emphasize here that the leading term in 
g'^{Q) for a given power of 1/Q^ is controlled by the relevant anomalous dimensions. 

We proceed by making an analogous asymptotic expansion of Ap{u) valid for u"^ G 
(A^, 00) with suitably large but otherwise arbitrary: 

Ap{u) = Appowcr{u) + Apcxp(?^) (10a) 

with 

Appowcr(«) - T. \:JL T.<^^n\T,p)[Hu')r, (10b) 
n=0 " iy=0 

where we have explicitly isolated in Apexp(w) all terms that vanish exponentially when 
— i> 00, such as the pole contributions to Ap{u) or terms containing the factor exp(— m/T). 
To simplify the notation we have chosen the units such that the running coupling has the 
form at one-loop level g'^{u) = 1/(6 In m^) (or equivalently, energy scales are measured in 
units of the relevant A-parameter) . 

This ansatz is sufficient to produce an asymptotic series of the form of Eq. (|^) to one- 
loop level. More generally, one could replace l/lnw^ with the full running coupling g'^{u) 
in the asymptotic sequence and generalize the method we are going to describe; the general 



strategy involved in this generalization can be found, for example, in Ref. |T^. Since we are 
presently only interested in one-loop calculations, we can regard g^{Q) and 1/lnQ^ to be 
proportional. It is easy to recognize that the sum over "n" is meant to match the sum over 
mass dimension in Eq. (P), while the sum over "z/" will match the sum over the order in 
perturbation series. The existence of anomalous dimensions in the OPE makes it necessary 
to introduce ^„ in the expansion for the spectral function. In the following we ignore the 
dependence of ^„ on "n" , as we have already ignored the fact that there generally exist more 
than one operator at a given dimension, and write ^ to avoid a too cumbersome notation. 
The complete notation will be restored when necessary. 

Before we proceed further, we wish to discuss whether it is possible that additional 
terms might appear in the expansion of Appoweriu) in Eq. ( p.Ob|) . In general, our asymptotic 
expansion procedure is powerful enough to exclude this possibility within the framework of 
the OPE. In fact, it allows to verify unambiguously that terms different from the ones already 
present generate, when substituted in the dispersion integral, terms that are missing in the 
OPE. Let us examine two specific examples that might be suspected to exist otherwise. First, 
dimensional arguments could suggest terms like g'^T/V^^- Our procedure shows that any 



6 



l/\/v? term in Ap{u) would generate a l/y^Q^ term in AK{iQ)^ which does not correspond 
to any known condensate, and it is therefore excluded. The second example is given by terms 
such as ln(T^/M^), which are naturally produced by elementary perturbative calculations of 
the spectral function: are such terms present in Eq. ( |10b| )? It is indeed true that such terms 
appear in the high-T expansion of the spectral function at fixed u"^. However, the expansion 
of the spectral function that is relevant for comparing to the OPE is a high-w^ expansion at 
fixed T. In the next section an explicit calculation in the Gross-Neveu model will illustrate 
the general fact that, contrary to the high-T expansion, the high-w^ expansion does not 
generate terms like ln(T^/M^). 

At this stage, we can already recognize a fundamental, and often overlooked, character- 
istic of the spectral function. If we insert the term Apoxp(^) in the dispersion integral, we 
only obtain pure powers of since the exponential convergence allows a naive expan- 

sion of the factor l/{u'^ + Q'^) (Watson's Lemma |TB|). Therefore logarithmic corrections, i.e. 
powers of g^{Q), come solely from the Appower(w) term. Because we know that the running 
coupling g^{Q) is always present in the OPE series, the term Appower(M) must be present 
in the subtracted SF, and it obviously dominates the asymptotic regime {v? oo). This 
fact immediately implies that only a finite number of moments of the subtracted SF can 
possibly be finite, i.e. the naive expansion of the factor 1/(m^ + Q^) is generally wrong, and 
that logarithmic corrections play a important role. 

A standard method to tackle the dispersion integral in the large-Q^ limit is the Mellin 



transform. However, the use of Mellin transform methods |]T5| is extremely cumbersome 
when logarithms appear in the denominator. Since inverse logarithms cannot be avoided in 
the spectral function, we need to resort to other means. 

We carry out the dispersion integral of the Appower(^) term by splitting the integral in 
Eq. (^ into three intervals: (0,A2), {X^,Q^) and ((5^oo). 

The integral over the first interval can be naively expanded in powers of since 
there are no convergence problems: 



AK,m ^ r du' ~ £ /'^ du'u^-Ap,^^M ■ (11) 

In particular, the leading term in l/Q^ is 

AK.itQ) ~ £ du" Appower(M) + 0{l/Q^) . (12) 

In the second interval, we use the asymptotic form of Appower given in Eq. ( |10b| ) and 
obtain (for instance by repeatedly integrating by parts) 

AK.m . r . f ,H|(4r(lng-)"-(lnA-)" 

, ^ (-1)' r(a, jn - I) In A^) - T{a, jn - I) \nQ') ] 

where a = 1 — u — ^, and T[a, z] is the incomplete Gamma function. Notice that when a = 



1 — ^ — ^ = the correct result for the term with Z = n is In y%- , which also corresponds to 
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its limiting value: limQ,^o[(ln Q^)" — (In A^)"]/a = In We then substitute in the second 
incomplete Gamma function that appears in Eq. ([T^) the following asymptotic expression 
which is valid for — > cxd and — 37r/2 < arg2; < 37r/2: 

-l)™r(l -a + m) 



r[a, z] = 



m=0 



z™r(l - a) 



The resulting expression for AK2{iQ) is 



1)' 



n = 



(,^ r(a,(n-/) InA^) ^(hiA^)" 



a 



(14) 



(15) 



-| 00 (i/) r CXD 00 

n2 n2n 2^ 2^ 



^ n,u=0 ^ 



-1) 



l+m 



r(m + l-a) (-l)"[lng 



2la 



i=0 m: 



fo (n - /)'"+i [In g2]m+i-a r(l - a) 



where again (lnA^)"/Q; and (lnQ^)"/a should be understood as In(lnA^) and In(lnQ^) when 
a = l-u-^ = 0. 

Let us notice that in the first line of Eq. (|T3|) any term with / such that the integral 
du^ u^'- Appoweriu) is finite can be formally identified as 



00 r 00 

EE 

l/=0'-n=0 



r(a, (n-/)lnA2 



{u) (In A 



2\a 



a 



/ (iM^M^'Appower(M) • (16) 



In particular, the leading l/Q"^ term in Eq. ([I5| ) is 



Ak2itQ) 



1 - (.jlXo^nlnA^ MilnX'r^ 



■)2 

00 



(17) 



r 00 00 



r(m + l-a) [ing 



I _L V n^"'' V V 



2la 



a 



We can similarly expand the integral over the third interval 



power 



U 



L M^. ^2^/ r(a,(n + 01ng^ 

\l+m 



-1 CO [u) 00 CO 
\ " " n \ " \ ^ 



r(m + 1 — a) 



± V n^'^) V V (-1)'+'" r(m + l-a) 



1_ 



(18) 



In the end we add the leading l/Q"^ contributions from Eqs. (p!2|), (|T^) and (|T^) to the 
corresponding contribution from the naive expansion of Apcxp and obtain 



AK{iQ) = / du 



Ap{u) 



+2E 



Ap [IngT^ - 

n2 + Q2^Q2+ Q2 

(-1)' r(2m + l + i/ + 



[lnQ2]. 



t [I ing^] 



2m+2 



+ 



' 1 \ 



(19) 
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where Ap is defined by 

pX^ poo 

Ap= du^ Appowcviu) + / du^ Apc^p{u) 

Since the integral /°° (ixx~^(lnx)~^ can be finite only if ^ > 1, we can use Eq. (|IBp to identify 
Ap as the zeroth moment of the subtracted spectral function, only when the asymptotic 
expansion of the subtracted SF in Eq. ( |l(Jb| ) yields ^ > 1. Otherwise, if ^ < 1, the zeroth 
moment is infinite, though Ap as defined by Eq. (pOD still exists. 



E. Exact Sum Rules 

We derive the sum rules by comparing the coefficient of l/Q"^ in Eq. (|T9]) 







and the corresponding coefficient in Eqs. 



Q-'Akm ~ cf)(K)A([Oi]«)[/(Q)]'' + O (^i^, [/(Q)]"+^j , (22) 

where we are only considering cases with n = 1 and d\ = 2. Furthermore, we are presently 
interested in exact sum rules that can be derived with one-loop calculations and, therefore, 
we only compare leading orders in l/lnQ^. Then there exist three possibilities, depending 
on the value of t], which can be calculated using the OPE. 

(1) If ?7 = 0, then ^ = 2 (or a bigger integer) and 

A^= r du^ Ap{u) = cf\K)A{[0,]^) . (23) 



This is the result one would get by naively expanding Q^/ {v? + Q"^) without worrying about 
the convergence of the moments and disregarding logarithmic corrections in the OPE. 

(2) If ?7 > 0, then ^ = 1 + 77 and 

Ap = / du^ Ap{u) = . (24) 
Jo 

(3) U 1] < (positive powers of InQ^ in Eq. (P)) or if the term InlnQ^ appears in Eq. (^, 
then ^ < 1 or ,^ = 1, respectively. This implies that 

roo 

Ap^ du^ Ap{u) = 00 . (25) 
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We remark that, even in this case when the moment is not finite, the asymptotic expansion 
is still well defined. It is nice to see that whether the zeroth moment of the subtracted SF 
exists is refiected directly through the leading power of g^{Q) in the OPE series. 

Our main results, Eq. (|23|) and Eq. (|^), can be expressed in physical terms as follows. 
The zeroth moment of a SF for a current J whose OPE series yields 1] > Ois independent of T 
and fi, while the same moment for a current with rj = changes with T and fi proportionally 
to the corresponding change(s) of the condensate(s) of the leading operator(s). Although 
cf\K) and A([(5i]k) can separately depend on k, their product must be independent of k, 
since the zeroth moment is independent of k. 



F. Discussion 



At this point several general comments are appropriate: 

(1) Our derivation relies on the fact that an asymptotically free theory allows a perturbative 
expansion at short distances, making practical the use of the OPE and of the RGE. We 
understand why only short distance physics is involved if we realize that the integral over 
frequencies reduces Eq. (|1]) to the ensemble average of the equal-time commutator of the 
currents. Therefore, results such as Eq. ( |23D and Eq. ( p4D are completely determined by the 
one-loop structure of the theory and the particular current under exam. 

(2) Flavor, or other non-dynamical quantum numbers, does not change the expansions at 
the one-loop level in an essential way. Therefore, one can derive analogous sum rules by 
using other kind of subtractions, instead of the one we adopted. One such example is given 
by the exact Weinberg sum rules (T = 0) in the chiral limit [l^ . 



(3) The derivation of sum rules for higher moments of the SF requires the complete cance- 
lation of all the lower dimensional operator terms, not just the leading g'^{Q). In particular, 
one also needs current quark mass corrections to the Wilson coefficients. Without appropri- 



ate subtractions, higher moments do not even converge [0. 

(4) It is essential to properly take into account the logarithmic corrections when deriving ex- 
act sum rules, since the logarithmic corrections not only dictate whether Ap satisfy Eq. (p3[) 



or Eq. (^) but they also control the very existence of Ap fl?]. This procedure is in sharp 
contrast with the usual QCD sum rule approach, where the convergence issue is by-passed 
by applying the Borel improvement by explicitly introducing a cut-off parameter (the Borel 
mass) . 

(5) We believe that the (T, /i)-dependent part of the leading condensate appearing in Eq. (^31) 
does not suffer from the infrared renormalon ambiguity. In fact, only the perturbative 
term Kq can generate contributions to the leading condensate that are dependent on the 
prescription used to regularize these renormalons. But Ko is independent of T and p: any 
prescription dependence cancels out when we make the subtraction in Eq. (|^). On the 
contrary, unless we generalize Eq. (0) and make other subtractions, sum rules that involve 
non-leading condensates are, in principle, ambiguous. 

(6) It is well-known that conserved operators are not renormalized and, barring anomalous 
violations, verify the same "classical" identities that can be derived at the tree level. This 
fact is also verified in the sum rules. In fact when both the currents and the operators 
are conserved [r] = 0) and we obtain the result of Eq. p3|), i.e. one can use the "naive 
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asymptotic expansion" to derive the sum rule. Conservation of the current alone is not 
enough to warrant a "classical" identity. 



III. GROSS-NEVEU MODEL 

In the preceding section we have derived sum rules valid for any asymptotically free 
theory. In this section we illustrate the procedure in the 1+1 dimensional Gross-Neveu 
model [0 in the infinite- limit. On one hand, we can derive the sum rules in Eqs. (p3D 
and (1^) by explicitly calculating Wilson coefficients, (3- and F-functions in the vector and 
pseudoscalar channels, following the general procedure discussed in section II. On the other 
hand, since this model is soluble, we can obtain the exact spectral function at arbitrary 
(T, /i), and then explicitly verify both that the sum rules are satisfied and that the asymptotic 
expansion of the spectral function has the form given in Eq. (|10b| ). Moreover, we can also 



calculate the condensates, and therefore explicitly check that the OPE really matches the 
asymptotic expansion of the exact spectral integral. 
The Lagrangian of the Gross-Neveu model is 

C = ^i^-dij + -^ [(^^)2 + , (26) 

where t/) is a two-component Dirac spinor and has A^-component in the internal space. 
Equivalently, we can write 

/: = ^(27.5-(T-^757r)^-^(a2 + 7r2), (27) 

where a and vr are auxiliary fields. The coupling constant g"^ is independent of and 
held fixed in the limit oo. This model is asymptotically free when D = 1 + 1, and 

the chiral symmetry is dynamically broken at T = and = to the leading order in 
1/A^ |18|. In the following we give the exact solution at A^ — cxd with finite (T, /i). Although 
formulas are explicitly considered in the symmetry breaking phase, they are also valid in 
the symmetric phase provided that the vanishing limit of certain condensates (such as the 
dynamical fermion mass) is properly taken. 

We wish to remark on the well-known fact that the limit N ^ oo here should not be 
interpreted as a starting point for an expansion of the model at finite A^, but rather as 
different model in itself, which is in fact the model we decided to use for the purpose of 
illustration. Moreover, there exist arguments |1T9| , ^ suggesting that this model in the limit 
A^ — > oo is actually more relevant to 3 + 1 phenomenology than the model with finite A^. 



A. Exact solution in the large- limit 

Since the procedure to obtain the exact solution in the large-A^ limit is rather stan- 
dard and some of the intermediate steps can be found in the literature, see for instance 
Refs. |[l8|j2l|JT9| j6|,^ , we only give the definitions and final results. 
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1. Gap equation and phase diagram 



The gap equation and the phase diagram can be derived from the first derivative of the 
effective potential 



(28) 

Lj„={2n+l)nT 



We introduce a momentum cutoff A and add the counterterm ^ln(4A^/K^), where k is the 
subtraction point. Then the renormahzed effective potential becomes 



^1 ^ In ^ + f(Vl±ZlIll^\ (29) 



where we have defined 



fix; T, /i) _ ^ ^ + g(^_^)/T ^ ^ ' (30) 

and, in the second equality, we have eliminated the subtraction point by introducing the 
dynamical fermion mass at (T, /i) = 0, i.e. M = Kexp(— vr/^f^). The effective potential is 
then obtained by integrating Eq. ppj): 

(1 + e~(^+'')/'^)(l + e"(^"'')/'^) 

(1 + g-{vW^+M)/T)(l _^ e-(v^F+^-/.)/T) 



(31) 

The gap equation is simply given by 



da 



= . (32) 

a=m{T,fi) 



The system possesses "metastable states" at those values of (T, fi) for which the gap equation 
has more than one minimum solutions: the solution with lowest y(m^(T, fi)) selects the true 
ground state. When (T, /i)'s are such that m(T, /x) = is the only minimum solution to the 
gap equation, these points in the (T, /i)-plane defines a second order phase transition line, 
which separate the symmetry-broken phase, i.e. m(T, /i) ^ 0, from the symmetric phase, 
i.e. m{T,fi) = 0. Taking the limit a = in Eq. (|29|), after some suitable rewriting, this 
critical line obeys the following equation 

1 ^9^-1 f"^ smh{k/T) r^ffi, T ) Ci-i^ 

2Jo k cosh[{k + fi)/2T]cosh[{k-fi)/2T] A/A;^^ ' '^^' ^ ^ 

whose numerical solutions are plotted in Fig. |l| in full line. The full line turns into a dotted 
line at the heavy dot, where another minimum (dash line) appears and becomes lower than 
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the minimum given by Eq. (|33|). Since the mass gap is finite on the dash fine (given by 
the solution of V^(cr^) = V^'(cr^) = with o" 7^ 0, where prime denote a derivative with 
respect to cr^), the dash fine is a first order phase boundary. The critical temperature at 
= is given by = {e""^ /'k)M k, 0.566933M. The "tricritical point" (the heavy dot in 
Fig. □) can be found by imposing that (a^ ^ 0) = V"{a'^ ^ 0) = 0: T3 ^ 0.318328M, /ig ^ 



0.608221M |19|. When T = the chemical potential at which the first order phase transition 



takes place is /ic = M/ \/2 . 



2. Bubble graph and correlator in space-like region 

The basic ingredient to build up the correlator is the fermion bubble graph. In the 
Euclidean region the pseudoscalar bubble, i.e. the free bubble graph between two J5 = tl>i'j5ip 
currents, can be conveniently expressed as 

n(a;) = {00) + AU{uj) , (34a) 

where 



' {uj)\^=iQ = ln(— ) + - / /^.> , 7^) (34b) 



TT m 

and 



AU{lo)1=,q = - dkf{E,;T,fi)^^^^. (34c) 

Here and in the following k, is always the renormalization subtraction point, m = m{T,fi) 
is the dynamically generated fermion mass, and Ek = \rrn?+W. We shall also use the 
short-hand notation 

1 POO fjU 

m^TT -^mf{E,-,T,t^). (35) 



For instance, the use of this notation makes possible to write the gap equation as 

dk 
Jo Ek 



In ^f{E,-T,f^)=n{{En. (36) 

m Jo Ek 



In the large- limit, the correlator in pseudoscalar channel is given by the geometric 
sum of the bubble graph: 

{JroJ.hA^) = n(^) (1 + j^^i^)) ■ (37) 

In the Euclidean region, we can then expand this correlator in the — 00 limit as 

13 



where we have introduced the running couphng constant 



l + (W2^^)ln(QV«:2)' ^^^^ 



and used the gap equation Eq. (|36|). 

We can recognize the first term in Eq. (^) as the perturbative contribution at (T, fi) = 0, 
{J5J5)^^\iQ), which corresponds to the term Ko{iQ) in the generic OPE expansion shown 
in section The other two terms in Eq. ( ^Sj) can be interpreted as contributions from 
the condensates. As we will verify later, in fact, m^/g'^ = ([('?/'V')^])t,/x5 while {{El)) — 
^m'^{{E^)) = {[9oo])t,ij, is the condensate of the energy-momentum tensor. We shall see 
that the perturbative contribution and the coefficients of the condensates are indeed the 
ones obtained by the OPE of J5J5, and therefore the expansion of the correlator shown in 
Eq. ( p8|) is of the form of the OPE. 

Note that in defining the condensates we have absorbed factors of in order to simplify 
the notation. 



3. Spectral function 



To obtain the spectral function we need the bubble graph in the time-like region cu^ > 

— Im n = eM X jtanhf ' ^' ^ ) + tanhf ) \ 40 

N ^ ^ ^ U^l-4myu^ I V 2T ^ V 2T ^ ' 

and 



where V J stands for the principle value of the integral. 
The complete spectral function has the form 

piu) = n6{uj' - ml)^^^^ + e{u^ - 4m2)p,„,t(^) • (42) 

The pole contribution comes from the bound state (the pion), whose mass is the solution to 
the equation 
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(43) 



in the region (0, 2m). This equation has in fact solution m.,^ = 0, and the coupling constant 
of this massless pion to its constituents is given by 
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The continuum part of the spectral function is related to the bubble Il{uj) in the time-like 
region through 

, ^ Imn(cj) 

~ (1 + g^ReU{uj)/Ny + {gHmU{uj)/Ny ' ^ ^ 

In order to get an idea of the dependence of the spectral function on (T, /i), we plot p{lj) 
for several typical values of (T, /i) , which are indicated in the phase diagram in Fig. I. In 
Fig. we show p{uj) for values of (T, /i) that go from the symmetry-broken phase to the 
symmetric phase through the second order line, while in Fig. ^ we show p{uj) as (T, p) go 
through the first order boundary. We immediately notice that the behavior of the spectral 
function changes drastically near the second order phase boundary. This is a clear reflection 
of the low energy critical phenomena related to the continuous phase transition. In fact, the 
pseudoscalar correlator (in the symmetric phase only, since it is then degenerate with the 
scalar correlator) at = can be regarded as the chiral susceptibility, and should diverge 
as \T — Tc\~'^ (with 7 = 1 in this model). At the boundary of the first order phase transition 
(Fig. ^), but still relatively close to the "tricritical point", the same qualitative deformation 
is present, since a weak first order transition shares certain qualitative features of a second 
order transition, even if we do not expect any real divergence. 

It is this kind of dramatic behavior of the spectral function near phase transitions that 
makes practical parameterizations nearly impossible. The chiral restoration in QCD is 
expected to be either continuous, or at least a smooth crossover, and the spectral functions 
in scalar and pseudoscalar channels should show a behavior qualitatively similar to the 
Gross-Neveu model, i.e. strong peak right above the origin near the phase transition. We 



shall comment more on this issue after we derive the sum rules for QCD in section IV 



4- Recovering the OPE series 

It is instructive to expand p{uj) in the limit uj"^ —>■ 00. We shall see that this asymptotic 
expansion is of the form assumed in section 0. Furthermore, once substituted in the disper- 
sion integral, it generates the OPE series of Eq. (^) according to the procedure developed 
in section |T[ 

To the leading order in we find 

|!lm(n(..)) ~ gMu){\ + ^ + 0(-l, e-H/^^)} (46) 

and 

1 + f^Re(n(.)) - + ^[2-^(1 - ^) - V(^)(M)rJ + 0{^^)] , (47) 



which leads to the following asymptotic form for the continuum part of the spectral function 

iV^---^") ~ 21+,V)/4 " 

+ ^(^)'a^oo])T4l + 0(/(Q))}. (48) 
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Here we only concern with the (?^(Q)-dependent terms and postpone the pure 1/Q^ terms 
(related to Apexp) to a later subsection. We recognize again the first term in Eq. as the 
asymptotic perturbative spectral function. Moreover, the (T, /i) -dependent part of Eq. 
has indeed the general form of Eq. (|10b|) . 



In connection with the comments we made after Eq. ( |10b| ), we point out that Eq. ( ^Hf ) 
has been obtained by expanding the spectral function in the limit ^ oo at fixed T. Had 
we made instead a high-T expansion (T —>■ oo) at fixed Q^, we would have obtained, for 
instance, ln(T^/co'^) from the third term of Eq. (|4l|). 

Upon substituting Eq. (H^) in the dispersion integral, and using the identity 



where z/ > 1 and we have kept only the 5'^(Q)-dependent part, we find the following series 
to the leading order in g'^{Q), 

y^Tw^^^ — N — ^'^^^ — — ^^^^ ~Q^^"' 

This result is in agreement with the OPE series of Eq. (|38D , which has been directly expanded 
from the correlator in the Euclidean region. 

The fact that no pure term appears in Eq. (|49|) implies that the zeroth moment 

of the subtracted spectral function vanishes: Ap = 0. The zeroth moment of the spectral 
function without the subtraction is not convergent, because of the contribution from the 
perturbative term, even if one can define it through a proper analytic continuation that 
yields p = 0. 



B. Direct solution in terms of the OPE 

We have just verified that in the Gross-Neveu model the exact correlation function can 
indeed be expanded in an asymptotic series precisely in the form of the OPE. Now we shall 
calculate the Wilson coefficients, /^-function and appropriate anomalous dimensions for the 
OPE, and verify that they match the coefficients and exponents of this asymptotic series. 
Furthermore, though it is beyond the scope of the OPE and RGE, we shall also calculate 
the condensates. In this way, we explicitly verify that the asymptotic expansion of the exact 
correlation function is identical to the OPE. 



1. One-loop calculation 

The one-loop /3-function of the 1+1 dimensional Gross-Neveu model has already been 
calculated in the original paper |T^: (3 = —g'^/2TT, i.e. according to our notation b = Xjl-n. 

The anomalous dimensions for pseudoscalar current J5 = ipi'^^ip and four-quark operator 
O4 = {ipipy can be calculated using the Feynman diagrams shown in Fig. ^ and Fig. ^b, 
respectively. According to the notation of Eq. (j^) we find 7jg = I/tt, and 704 = 2/71. The 
energy-momentum conservation makes the anomalous dimension for ^00 vanish, i.e. ^0^^ = 0. 
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The relevant Wilson coefficients corresponding to operators O4 and Oqo can be calculated, 
to leading order in g'^{Q), using the Feynman diagrams shown in Fig. |]a and Fig. ^b, 
respectively. The resulting coefficient are C04 = —2g'^{Q) (z/q] = 1 according to notation in 
Eq. (D) andc,,„=4(j.(°; = 0). 

These results lead to the exponents 7704 = 1 + (27J5 — 704)726 = 1 and rjg^g = (27 — 
7eQQ)/26 = 2, which are precisely what we found in Eq. (^Sp. 

These explicit one-loop calculations exactly match the solution in Eq. p9|). 

At last a comment on the one-loop calculation: this result is exact to the leading order in 
1/N and only involves insertion of the quark-bubble chain. Since the quark-bubble behaves 
like oc InQ^, this insertion is equivalent to substituting the coupling constant with the 
running coupling constant g'^{Q)- So one can easily identify those sets of Feynman graphs 
whose sum leads to the solution of the RGE in Eq. (H). 



2. Condensates 

The OPE itself does not specify how the relevant condensates are calculated, since the 
OPE is only a RGE-improved perturbative procedure. The calculation of the condensates 
can be done only in a non-perturbative context. Here, we carry out this calculation using 
the expansion. 

The bare quark condensate {ipip)T,ii is given by the dynamical quark tadpole graph 

mr., - i-^T) f ^ E IV^ ^ 2»r g E . (50) 

Using the standard contour integral technique to carry out the Matsubara frequency sum, 
we obtain 

(^^)T,, = m((E°))--ln(-). (51) 

The renormalized quark condensate is obtained by replacing 2A with the subtraction point 
K in the above equation 

{mr,, = m{{El)) - ^ In(^) = . (52) 



In the last step we have used the gap equation Eq. (p6D. Notice that the renormalization of 
this condensate is not independent of the one carried out for the effective potential. In fact, 
the cr-field is just an auxiliary field: a = —g'^tpip. Therefore, the gap equation gives us not 
only the expectation value of the cr-field, but also of ipip: {cr)T,fi = ^iT,fi) = —g^{[ilJip])T,^i- 
So the renormalizations of ipip and g^ are related. 

The four-quark condensate factorizes in the large- limit, and we find {[{'ipipy])T n = 

The bare kinetic energy expectation value is given by 



^ dk ^ i dk ^ {iun + fi) 



2 



7 — m J 2n ^ {iujn + /x)^ — {k"^ + rn?) 

(53) 



17 



After removing the (T, /i)-independent volume and quadratic divergences, and introduced a 
subtraction point for the logarithmic divergence the renormalized kinetic energy becomes 

2 

([^^7o9o^])T,, = {{El)) - — In(-) . (54) 

After subtracting the trace term from the kinetic energy, we finally obtain the expectation 
value of the traceless energy-momentum tensor 

([^oo])t,, = ([^no9o^ - -m^^])T,, = {{El)) - —{{El)) . (55) 

In the symmetric phase (m = 0), the expectation value of the energy- momentum tensor has 
the very simple form at = 

([M)t,,=o = (56) 

One can verify that inserting these condensates in the OPE, whose coefficients and 
exponents we have already calculated, reproduces the asymptotic expansion of the exact 
correlator in Eq. (pSl) . 



C. Explicit verification of the exact sum rule 

Let us first notice that the sum rule / duj"^ /\p{uj) = converges very slowly. In fact, the 
subtracted spectral function has the leading behavior 

Apiuj) ^ Ap,{uj) + 0{-) , (57) 



where 

^P'^"^^ - -^^['7(^4. + TRTi / • ^^^^ 

Since g'^{uj) = n/ln^u/M), the leading behavior in Eq. (|58|) is Ap2{uj) ~ [c<j^(lnu;)^]~\ and 
it is clear that the sum rule only converges logarithmically, i.e. lim^2^oo J^^ doo'^Ap{uj) ~ 
l/lnv42. 

Moreover, it is only thanks to the the logarithmic corrections that the sum rule is finite. 
Had we disregarded this corrections, we would have found Ap2{uj) ~ uj'"^, and the zeroth 
moment of the subtracted spectral function would not exist. 

Since the exact spectral function cannot be integrated analytically, we verify the exact 
sum rule numerically. We take advantage of the fact that the integral of Ap2{uj) in Eq. (|58D 
can be done analytically to improve the numerical convergence of the integral. We rewrite 
the sum rule as 

rA^ roo roo 1 

duj^ Ap{iu) = - j^^ duj^ Ap{Lo) = - j^^ duj^ Ap2{uj) + 0{—) , (59) 
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where the integral of Ap2{uj) can be done analytically and we obtain 



I d.^ApM = -{2A,n'(-taT- — ^l) + } + O(^) . (60) 

This formula can be easily verified numerically with high accuracy, since the error is now 
rather than 1/ In A. It is reminded that Ap in the above equation should be obtained 
from the complete SF in Eq. (^2]), including both Appower and Apexp- 



D. Correlator in vector channel and baryon number susceptibility 



There are several reasons to study the correlator in the vector channel in the Gross- 
Neveu model. First, the vector channel gives us an example where the zeroth moment of 
the subtracted SF does not exist, but the asymptotic procedure can still be carried out. 
Moreover, in connection with the fact that the vector current is conserved, this spectral 
function is not affected by logarithmic corrections, and we can use it to illustrate the use of 
the Mellin transform method to obtain the asymptotic expansion. In addition, very much 
like in the QCD case, we shall see that the baryon number susceptibility changes drastically 
near the chiral restoration transition, even if we do not expect any real singularity associated 
with critical phenomena in the vector channel. For convenience, we only consider the case 
of p = 0. 



1. Vector correlator in Euclidean region 



The Feynman graphs that contribute to the vector correlator are depicted in Fig. ^. Note 
that if there were no mixing between vector and pseudoscalar channel, which in fact is only 
present in 1+1 dimensions, the sole contribution would come from the first graph. The total 
contribution from all the graphs in Fig. is 



{JaJ )tA'>'Q) = / = tanh W ^„ . „ + 2m 



^ ^ ^Q^iQ^ + Am^) V + 4m2 Jo Ek El + gV4 J ' 

(61) 

We define the spectral function in vector channel as 

(/„/"k.w)=i{i-/;<i.^^}. (62) 

Then we find the following explicit expression for p 

P(^) = / .. . . . (l - eM/2T + i ) ■ (63) 



It is clear that leading asymptotic behavior of the spectral function is p{u) ~ 1/t^^, which 
implies that the zeroth moment of the spectral function does not exist (it is divergent). In 
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the section H we have shown that, corresponding to this divergence, we must find negative 
powers of g'^{Q) (positive power of \n{Q'^)) in the OPE series. 

When there are no exphcit logarithms in p{uj), it is easier to make the asymptotic ex- 
pansion by using the Mellin transform method, and in particular the convolution property 



of the Mellin transform ||15|| , 

roo 1 rc+ioo 

/ dtf{t)h{t) = — dzM[h;z]M[f; 

Jo Zm Jc-ioo 

where M[h; z] is the Mellin transform of h{t), defined by 



l-z] 



M[h- z] 



dte~^h{t) . 



depends on the asymptotic behaviors of h{t) at t 



The value of c in Eq.( 
at t ^ respectively; details can be found, for instance, in Chapter 4 of Ref. |T5 



(64) 
(65) 

oo and f{t) 



Now let us illustrate how to use this method to expand the power part of the spectral 
function (the part containing exponential dumping factors can be expanded naively) 

Am'^6{uj'^ - Arn^) 

Ppower K.^ ) I ! 

Juo'^{uj'^ — 4:1V?) 

whose Mellin transform is 

M[ppower; z] = {Am'y^Vil - z)/V{?>/2 - z) . 
Using the convolution, Eq. (|64D, and the Mellin transform of 1/('U^ + Q"^) 



(66) 



(67) 



M[(m2 + Q 



2\-l 



;1 



7r(g')-7sin7r(l-2) 



(68) 



we find (in this case c in Eq.(|64D obeys < c < 1 



du' 



Ppov 



Res 

n=l 

4m^ 



4m2Nz0Fr(l 



IT 



g2 J r(3/2-z) sin7r(l-2) 



ln( 



Q' 



l + ln(«j)]+0( 



(69) 



Similarly to the pseudoscalar correlator, it is possible to check this result in two ways. We 
can directly expand the power part (temperature independent part) of Eq. ([^), and we can 
derive the OPE {'jv = and C04 = 8g'^{Q) for the l/Q^ term). 

As promised, there indeed appears a positive power of ln((5^) in the OPE series. At the 
same time, the singularities of M[ppower; z] at z = positive integer hints the fact that all the 
non- negative integer moments of Ppower do not exist. On the other hand, the asymptotic 
expansion in Eq. ( |69D is well defined. Even though the vector current is conserved, due to the 
fact that the anomalous dimension of O4 is not zero, the correct result can not be obtained 
by a naive asymptotic expansion in this case, in accordance with the remark number (6) at 
the end of section II. 

This same result could have been obtained also by the general method developed in 
section ||, but it should be clear by now that the use of the Mellin transform, when possible, 
is more straightforward and require less labour. 
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2. Baryon number susceptibility 



The baryon number susceptibility (xbn) is defined as the vector correlator in the limit 
g2 ^ 0. From Eq. (|61|) we find 

XBN = vr(J.J")HzQ^0) = -l + 2m2y^ ______ (70) 

Notice that since we have contracted the Lorentz indices when defining our vector correlator, 
Xbn is nonpositive definite. It is also possible that our definition of xbn differs from others by 
an additive constant, which is irrelevant to the temperature dependence we are considering. 
In the symmetric phase, the limit m — > should be take with some care: the result is 
Xbn = 1, when T > Tc. The complete temperature dependence of xbn is displayed in Fig. ^b. 
As expected, there is no divergence or critical phenomena in this channel. Nevertheless, we 
do see a sudden rising of xbn in correspondence to the relatively rapid drop of the dynamical 
quark mass, which we show in Fig. |^a for comparison, near the phase transition region. 

IV. APPLICATION TO QCD 

After a lengthy detour to the Gross-Neveu model let us now back to the main interest 
of our work, to derive exact sum rules and to explore the corresponding phenomenological 
consequences in QCD. As we will see, QCD shares many of the qualitative features observed 
in the Gross-Neveu model, implying the simple fact that any drastic change induced by 
a phase transition or rapid crossover would necessarily be reflected through the spectral 
function in proper channels at low frequencies, independent of the details of the model. 



A. Sum Rules 

Let us consider four correlators of mesonic currents in QCD. The correlation between 
scalar currents and the one between pseudoscalar currents involve non-conserved operators, 
Js = ipip and Jp = ^"^^ip, whose anomalous dimensions are = = l/47r^. One the 
other end, correlators between vector currents and axial-vector currents involve conserved 
quantities, Jy = 'ip'~i^ip and J a = ^T/^Ts'?/', whose anomalous dimensions vanish (7^ = 7a = 
0). 

In QCD there exist several dimension-four operators, but all of them have non-positive 
anomalous dimensions, therefore the two correlators between nonconserved currents have 
Vn > (27j — 7n)/2& > and Eq. (p^) applies, i.e. the zeroth moments of their SF's are 
independent of T and /x. 

On the other hand, the two conserved currents have, in correspondence with conserved 
operators = and a generalization of Eq. applies. In this case in fact there are 
three dimension-four operators with zero anomalous dimension. Two of these operators are 
Lorentz scalars: mipip and agG"^ = astT{G^iyG'^'^), while the third is the energy-momentum 
tensor: 

1 - ^ ^ 1 

Ofiu = -ij{i^iD^ + '^u^D^ - mg^^)i) + 2tr(G'^„G'/) - -g^.yG'^ . (71) 
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The sum rules for the vector and axial-vector currents (with the Lorentz indices contracted) 
are: 

/ du'Apiu) = aA{[mi,^p])+ + 8A([goo]) , (72) 

with a = 6 and a = — 10 for vector and axial- vector cases, respectively. 

Obviously, the general remarks at the end of section II apply to the sum rules derived 
above. However, there is a new feature of QCD which may not be shared by all asymp- 
totically free theories. It is known that QCD OPE series has explicit instanton induced 
corrections that cannot be related to condensates of some local operators. But these exact 
sum rules are not explicitly affected by this problem, although the value of the condensates 
certainly have instanton contributions. The reason is that the instanton singularities in the 
Borel-plane are located on the positive axis starting at Svr^, which in turn implies that the 
explicit contribution of instantons to correlation functions are of higher order than l/Q^ ||2^ 
in the OPE series. 



B. Phenomenological Consequences 

Finally, let us discuss some of the phenomenological consequences of these exact sum 
rules. The power of these exact sum rules rests on their generality, in the sense that there 
is no assumptions on the functional form of the SF's and they can be applied anywhere in 
the (T, /x)-plane. 

First let us consider the pseudoscalar channel. The sum rule Ap = implies that, in 
the broken-chiral-symmetry phase, the change of the pion pole induced by T or /i is exactly 
compensated by a corresponding change of the continuum part of the SF. 

Next let us consider the scalar correlation function at zero frequency, i.e. the chiral 
susceptibility (a measure of the fluctuation of the chiral order parameter), 

x{T,p) ^ fd'xeixo){[Jsix),Jsm)T,, = rdu' (73) 

which diverges when (T, p) approaches the phase boundary, provided the chiral restoration is 
a continuous transition. The divergence of the chiral susceptibility near the phase transition 
can only be produced in Eq. ([73|) by singularities very close to the origin; singularities not 
near the origin are not compatible with the sum rule Ap = 0. Thus, when approaching the 
phase boundary in the (T, /i) plane, the threshold of the spectral function vanishes (since 
there is no massless pole in the chirally symmetric phase), and a strong peak develops right 
above threshold. In the chirally symmetric phase the pseudoscalar and scalar channels are 
degenerate, therefore we expect the same behavior for the pseudoscalar SF. This strong peak 
in the pseudoscalar and scalar SF's, which is intimately connected with critical phenomena of 
a diverging susceptibility and correlation length near the phase transition, can be interpreted 
as some kind of quasi-particle, thus confirming the qualitative picture, originally proposed 
in the context of the Nambu-Jona-Lasinio model [^|, of the appearance of soft modes near 
the chiral phase transition. 
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We wish to remark that the sum rule has been used only to exclude the logical possibility 
of the SF developing singularities at finite energies. But the sum rule is not a necessary 
condition for the appearance of a peak. The presence of a second order phase transition 
(infinite correlation length) is the real physical cause of the peak. Nevertheless, one should 
not assume that if the phase transition is not second order the zero-T or perturbative spectral 
function is a good approximation to the real spectral function. Even if the chiral restoration 
turns out to be a cross-over or weak first order transition (finite but large correlation), as 
the lattice data seem to indicate ||2^, we still expect the same qualitative features, though 
less pronounced: a peak develops but it does not actually diverge at the transition. Our 
calculation in the Gross-Neveu model is a clear illustration of this situation. While Fig. (0) 
shows the divergent peak of the spectral function near the second order boundary, we see in 
Fig. (^ that the spectral function is also strongly peaked near the first order transition as 
long as we are not too far away from the "tricritical point" . 

We believe that similar results hold in the vector and axial-vector channels, even if our 
conclusions cannot be as strong as in the previous case, because of two main differences. The 
first is that now Ap 7^ 0. This problem is not very serious, since we only need that Ap is not 
singular crossing the phase boundary. And this can still be argued by means of the sum rule 
in Eq. ( [7^ ) and the fact that, for a continuous phase transition, the changes of the thermal 
energy A([6'oo]), and of both condensates A([m'?/'?/']) and A([asG^]) should behave smoothly 
across the critical line. The most serious difference is that, in these channels, we cannot 
argue on physical grounds that the corresponding susceptibility diverges. Nonetheless, there 
exist lattice simulations p6| showing that the so-called quark number susceptibility rapidly 
increases in the transition. Thus it is plausible that also in the vector channel the SF has 
a vanishing threshold, and accumulates strength just above it. A similar argument in the 
context of the Nambu- Jona-Lasinio model can be found in Ref. 1271 . 



V. SUMMARY AND OUTLOOK 

We used the operator product expansion and the renormalization group equation to 
derive exact sum rules at finite T and fi valid for asymptotically free theories. Our derivation 
explicitly show that logarithmic corrections cannot be neglected. We found that, depending 
on the theory and on the current under study, the zeroth moment of a spectral function is 
either independent of T and fi, or its change is related to the corresponding changes of the 
condensates of operators of lowest dimension. 

In particular the zeroth moment of the scalar and pseudoscalar mesonic currents in QCD 
are independent of T and /i. As a consequence any change in the strength of the pion pole 
must be exactly compensated by a change in the continuum contribution. We also infer that 
spectral functions in the scalar and pseudoscalar channel should dramatically change near 
phase transition. We also find that the zeroth moment of the vector and axial-vector mesonic 
currents in QCD changes with T and /i, and these changes are related to the corresponding 
changes of the condensates. 

Due to their generality, these exact sum rules strongly constrain the qualitative shape 
of SF's, in particular, near phase transitions. Since it appears that a strong deformation of 
the spectral function from its counter part at (T, /i) = and perturbative cases is the most 
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likely scenario near phase transitions, the information carried by these exact sum rules is 
particularly welcome. We urge whoever parameterizes a spectral function, e.g. in QCD sum 
rule type of calculations or to interpret lattice simulation data, to incorporate these exact 
constraints. 

We have also illustrated in great detail the derivation of the sum rules, and confirmed 
their validity, in the Gross-Neveu model, where we can also calculate the exact spectral 
function. Even if this model calculation is meant mostly as an illustration of many delicate 
issues of the derivation (such as convergence, asymptotic expansion, Mellin transform, high 
momentum expansion vs. high temperature expansion, etc.), it is nevertheless comforting 
to find that all our general expectations about the spectral function are explicitly verified 
in the model. 

In the future, we plan to generalize these results to baryonic currents and to analyze the 
phenomenological consequences of these sum rules in greater detail. 

This work is supported in part by funds provided by the U.S. Department of Energy 
(D.O.E.) under contracts #DE-AC02-76ER03069 and #DE-FG06-88ER40427. 
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FIGURES 



FIG. 1. Phase diagram of the Gross-Neveu model in 1+1 dimension. The continuum line 
represent the second order phase boundary in the (T/Tc, /x/M)-plane, while the dashed line is the 
first order phase boundary. The "tricritical point" , where first and second order transitions merge, 
is shown as a heavy dot. The dotted line separates the values of (T, n) for which the symmetry is 
always broken from those values of (T, ^) for which a "metastable symmetric phase" is possible. 
Dots show those values of (T, /i) at which we plot the spectral function in Fig. ^ and Fig. ^. Their 
coordinates are: a = (0.352, 0.4), b = (0.848, 0.4), c = (0.875, 0.4), d = (1.763, 0.4), e = (1.175, 0.4), 
/ = (1.763, 0.4), g = (0.353, 0.63), h = (0.353, 0.65) and k = (0.353, 0.85). 

FIG. 2. Spectral function of the 1+1 dimensional Gross-Neveu model in pseudoscalar channel. 
We show several values of (T, /i) that go from the broken phase to the symmetric phase passing 
through the second order phase transition. The labels a, b, c, d, e and f correspond to the labels 
on the points in the (T/Tc, ;u/M)-plane shown in Fig. |l[ The upper panel shows spectral functions 
in the broken phase. The pion poles are indicated by arrows, and are slightly displaced from their 
true position at u; = for vision convenience. The length of arrows represents the magnitude of the 
corresponding residues. The lower panel shows spectral functions in the symmetric phase, where 
there is no pion pole. 

FIG. 3. The notations are the same as in Fig. ||, except that now {T,fi) goes through the first 
order phase boundary. The coordinates of the labels a, g, h and k are given in Fig. |^. 

FIG. 4. Feynman diagrams for calculating the anomalous dimensions of (a) the composite 
operator J5, i.e. 7jg, (b) and of the composite operator O4, i.e. 70^. 

FIG. 5. Feynman diagrams for calculating Wilson coefficients: (a) C04 and (b) cggp. 

FIG. 6. Feynman diagrams contributing to the vector correlator. 

FIG. 7. (a) the dynamical quark mass and (b) the baryon number susceptibility as functions 
of temperature. 
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